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We study the photocurrent response of topological insulator surface states to circularly polarized 
light for arbitrary oblique incidence. We describe the surface states within a Dirac model, including 
several perturbations such as hexagonal warping, nonlinear corrections to the mode velocity, and 
applied magnetic fields. We find that the photogalvanic current is strongly suppressed for the usual 
orbital coupling, prompting us to include the weaker Zeeman coupling. We find that the helicity- 
independent photocurrent dominates over the helicity-dependent contributions. 



Introduction. — Topological insulators (TI) recently 
emerged as a central theme in condensed matter physics.^ 
The intense interest in this new state of matter is rooted 
in their unique properties .'^HH In addition to a band 
gap, TIs have conducting surface states with remark- 
able properties. They are protected against backscat- 
tering by time-reversal symmetry, and are helical: each 
surface-momentum state possesses a unique spin direc- 
tion. The unique properties of the surface are re- 
sponsible for their exotic electromagnetic properties,'^ 
and might be used to realize topological superconduct- 
ing phases hosting Majorana modes when brought into 
contact with s-wave superconductors.-'^'^ Bismuth-based 
compounds were among the first materials predicted to 
be three-dimensional TIs,'^^^^^^ a prediction verif ied ex- 
perimentally by angle-resolved photocmissiorPISHIII and 
scanning tunneling spectroscopy.^^^^ISl 

The surface states also exhibit exotic optical prop- 
erties. Gapped surface states are predicted to cause 
giant Kerr and Faraday rotations of polarized light .^221 
The helical nature of surface states is expected to make 
their photocu rrent response to electromagnetic radiation 
rather unique .f^il^Il As illustrated in Fig.[T[a), due to spin- 
selection rules, one might expect that circularly polarized 
light excites the surface states anisotropically around the 
Fermi surface, thus inducing a dc electric current. Such 
photocurrents constitute an interesting probe of surface 
states in TIs (see, e.g., Ref. [53]). An important moti- 
vation for investigating the photoresponse of TIs is that 
helicity-dependent currents are expected to emerge solely 
due to the surface, as they require the breaking of in- 
version symmetry precluding a second-order photogal- 
vanic effect.!^ Also, bulk electronic states in TIs consist 
of Kramers pairs, which rules out a helicity dependence 
of the photocurrent as welLl^^l 

Here we investigate TI surface photocurrents within 
a minimal model of the surface states motivated by 
Bi2Se3, and obtain rather surprising results. The inci- 
dent light couples to the surface electrons in two ways: 
through orbital minimal coupling (Peierls substitution), 
p — >■ p — eA, and through the Zeeman energy. Usu- 
ally, the orbital effect is expected to dominate strongly, 
and the Zeeman coupling is neglected. Here, however, 
we find that the orbital component of the photocurrent 



vanishes for the simplest model of a perfect Dirac cone, 
as all surface-electron spins lie in the plane of their mo- 
tion, and light excites carriers isotropically around the 
Dirac cone. The orbital coupling can induce photocur- 
rents only when perturbations of the ideal Dirac c one are 
included: hexagonal warping of the Fermi surface f^^'^HIIl 
an external magnetic field; and a mome ntum-dependent 
correction to the Fermi velocity.l^^l^lllll Since these per- 
turbations are quite small, we also include the Zeeman 
coupling to the incident light in our analysis. Surpris- 
ingly, the Zeeman coupling is responsible for the domi- 
nant contribution to the surface photocurrent response, 
which we find to be helicity independent, linear in the 
Zeeman coupling, and to flow against the direction of 
propagation of the light. The helicity-dependent pho- 
tocurrent, suggested by the simple mechanism illustrated 
in Fig. [ij^a), is found to be very small, i.e., quadratic in 
the Zeeman coupling. 

Model. — We perform our analysis within a minimal 
model of a TI surface. Let the surface lie in the xy- 
plane, with radiation incident at an arbitrary angle as 
illustrated in Fig. [T] We assume that the energy huj of 
the incident radiation is such that the excitation takes 
place solely within the Dirac cone located within the bulk 
band gap (i?gap = 0.3 eV for Bi2Se3^^ ^^). As usual, we 
neglect the momentum change in the optical transitions. 
Motivated by the (111) surface of Bi2Se3, we consider the 
model Hamiltonian 

H = v{pxay -pyOx) + ^(j3+ +P^-)(Tz - g^isBax, (1) 

which includes cubi c war ping 

X^smm a correction to 
the Fermi velocity A,'22l2ll and an external magnetic field 
B in the a:-direction. Here v = vp + Ap^ with the Fermi 
velocity vp, g denotes the g- factor, /is is the Bohr magne- 
ton, p± = px i: ipy = pe^"^, and a± = i: iayW^ This 
Hamiltonian is partice-hole symmetric with H\p, ±) = 
±E\p, ±). We consider circularly polarized light incident 
onto the sample at an azimuthal angle ip from the nega- 
tive y-axis and at a polar angle 9 from the positive z-axis 
(cf. Fig. [T]). For If ~ and left-circular polarization, the 
vector potential is given by 

A{t) = Ao[cos{kr — Lot)x — sin(fcr — Lot){x x k)]. (2) 
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mensions can be written as 



(4) 



where the sum is over positive-energy states, rip {np _^^) 
is the distribution function (equiUbrium distribution 
function) of momentum state p in the positive energy 
band, i'p,+ is the velocity of a particle in state p in the 
positive energy band, and the factor of 2 accounts for 
particle-hole symmetry. Assuming that momentum re- 
laxation occurs on a much faster timescale than energy 
relaxation, the steady-state solution of the Boltzmann 
equation in relaxation time approximation gives 



FIG. 1. (color online) (a) Basic idea for the generation of 
a helicity-dependent photocurrent. Dark red arrows indi- 
cate spin direction, filled (empty) circles indicate electrons 
(holes). Circularly polarized light induces spin-dependent 
transitions, exciting electrons asymmetrically in fc-space. (b) 
Illustration of the warping effect on the Fermi surface. Lighter 
colors/shades correspond to larger absolute values of energy. 
Dark red arrows indicate spins, circles with -|-(— ) indicate 
positive (negative) spin z-component. (c) Illustration of the 
orientation of the surface, the direction of the incident light 
and the resulting main contribution to the current. 



While we focus on the response to circularly polarized 
light, we find that the surface also exhibits a photocur- 
rent for linearly polarized radiation. For ip = 0, the di- 
rection of propagation k lies in the yz-p\a,ne at an angle 
6 from the positive z-axis, such that k — sin 9y ~ cos 6 z 
with 6le[0,7r/2] (see Figjl]). 

We include both orbital and Zeeman coupling of the 
light to the surface electrons in the coupling Hamiltonian 



where Tp is the momentum relaxation time and T is the 
transition rate from state |p, — ) in the lower band to state 
|p, in the upper band which can be calculated usinj 



Fermi's Golden Rule. For T = 0, (r7.° _ 



in Eq. 5 



is only nonzero if the chemical potential lies between the 
energies of the two states participating in the transition 
(cf. Fig. [l]). This gives a condition for the minimum 
photon energy required to induce transitions. 

Assuming huj/2 > for the photon energy, the cur- 
rent density becomes 



Aner, 



H' = —eVp ■ A — .9sM-b(V X A) • <t. 



(3) 



J2vp,+ \{P,+\H'_\P,-)\'S{2E - hw), (6) 



where we used that i?p,+ = —Ep^^ = E due to particle- 
hole symmetry. 

The calculation of the integrand in Eq. ^ is expedited 
by using projection operators P± onto the two bands, 



where m is the electron's mass, a the vector of Pauli ^p,+ I(Pi )P — tr 



matrices, and Vp = OH /dp. H' encodes both photon 
absorption and emission processes, corresponding to cre- 
ation or recombination of an electron-hole pair, respec- 
tively. It is, therefore, useful to write H' = H'^ -f H'_ 
with H'^ (X A± = Aoe^'"' describing the emission and 
absorption of a photon. 

The relative magnitude of Zeeman and orbital cou- 
pling can be estimated as ~ hk/mvp = huj/mvpc. For 
a photon energy of fiw = 0.1 eV and a Fermi velocity of 
Vp = 5 ■ 10^ m/s,l^this ratio is of order ~ 10~^. Even 
though the orbital coupling clearly provides the domi- 
nant excitation process, we find that it produces no net 
photocurrent without an applied magnetic field. For this 
reason, we include the Zeeman coupling and find that 
the leading contribution to the current for realistic val- 
ues of the magnetic field is given by an interference term 
between Zeeman and orbital coupling. 

Photocurrents. — For a particle-hole symmetric system 
without spin degeneracy the current density in two di- 



P^H'P+H'_P- 



dH 
dp 



where we used that v. 



value of an operator can be written as 



(7) 

-Vp._, that the expectation 
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2E J 



(8) 



and that the transition |p, +) — > |p, — ) (|p, — ) — > \p, +)) 
happens via emission (absorption) of a photon and is 
therefore mediated by the coupling involving Aj^ (^-)- 

In addition to being careful with the definition of the 
velocity, we must also make sure that the argument of 
the (5-function in Eq. ([6| contains the perturbations of 
the pure Dirac spectrum. To expand the argument of the 

oo 

(5-function, we use 6{2E - huj) ^ 4:E f ^gM^E^-f^^'^^) , 

— oo 

After expanding in the parameters A, A, and B to first 
order in A and A and second order in B and performing 
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TABLE I. Helicity-dependent corrections to the current in- 
duced by various perturbations, vz = vz/vf ~ 10~^, 
A = A(MV(4) ~ 10"', A = A(MV(4) ~ 10-^'^and 
B = g^sB /(hbj) ~ 10~*-B/T are dimensionless parameters 
and C is given in the text. 
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TABLE IL Hehcity-independent corrections to the current 
induced by various perturbations. Parameters as in Tab. |l] 



the angular integral in Eq. we find (see App. [A| ) 



da 



dp p 



= (0) 



X e 



ia[4(pVFf-{hwf] 



(9) 



where ■* are functions of momentum and A, A, and B. 



The integral over a is simplified by writing a - 
for the factors of a in the brackets, and then first carrying 
out the a-integral. The remaining integrals can now be 
easily done, since the complicated angular dependence of 
the eigenstates has been eliminated. 

Results for ideal Dirac spectrum. — For the simplest 
model of a perfect Dirac cone without external fields, 
i.e., A = A = i? = 0, and oblique incidence, we find that 
no net charge current is induced by pure orbital coupling. 
Although the coupling between vector potential and elec- 
tron momentum leads to much larger excitation rates, the 



transitions take place isotropically around the Dirac cone 
and no net charge current is induced. However, including 
the small coupling between vector potential and electron 
spin, we find that currents are generated perpendicular 
and parallel to the plane of incidence with different po- 
larization dependencies. We find that the current is given 
by 



C 

= - —Vz sin 6 [y' - vzx') , 



(10) 



where C 



rotated coordinate system such that A is incident in the 
y'z-plane. j^y^ is helicity independent, results from an 
interference effect between orbital and Zeeman coupling, 
and can also be induced by light which is linearly polar- 
ized perpendicular to the plane of incidence (S-polarized) 
(see App.|b|. The smaller current component j^? is he- 
licity dependent, i.e., it changes sign for uj — > — w, results 
from pure Zeeman coupling, and will not be present for 
linearly polarized light. 



, and x' , y' define a 



The result given by Eq. ( 10 ) could have been antic- 
ipated by symmetry arguments. The (111) surface of 
Bi2Se3 h as, among others, a mirror axis along the y- 
direction.lSSl When the light is incident with (p — 0, the he- 
licity of the vector potential changes sign under this mir- 
ror transformation, the current in the a;-direction changes 
sign, and the current in the y-direction remains invariant. 
Thus, helicity-dependent currents are only allowed in the 
x-direction while helicity-independent currents must flow 
in the y-direction. Since the system is rotationally sym- 
metric for A = _B = 0, the only directionality is provided 
by the vector potential and the currents will rotate ac- 
cordingly. In addition, the interaction Hamiltonian for 
93 = is given by 

H'_ ^ vf {(^y — i COS 9ax) — vz {(Jx + i cos Ouy + i sin Oa^) ■ 

(11) 

For normal incidence {9 — 0), the orbital coupling is pro- 
portional to the spin raising operator in the z-direction, 
cr+ = (Tx + icTy, while for oblique incidence it involves 
a sum of spin raising and lowering operators in the z- 
direction because Cy — i cos 6ax ^ 0-+ ( 1 -I- cos 6) — a-{l — 
cos 9). Since for A = all spins lie in plane, electrons 
are excited isotropically around the Dirac cone and the 
orbital coupling by itself cannot generate a net current. 
For normal incidence the same argument even excludes 
currents induced by the Zeeman coupling. For oblique in- 
cidence, however, the Zeeman coupling involves, through 
ax + i sin 9a z, a sum of spin raising and lowering opera- 
tors in the y-direction. Since the lowering operator has 
the larger coefficient, the Zeeman coupling will preferably 
excite spins with momentum in the negative x-direction 
generating a current in the x-direction (cf . Fig. [T]) . Simi- 
larly, we find that in the interference term the spin lower- 
ing operator in the a;-direction a^+iay dominates, prefer- 
ably exciting electrons with momentum in the positive y- 
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FIG. 2. Current as a function of polar angle (fi of vector 
potential, hd (hi): helicity dependent (independent). 9 = 
0.98 and the other parameters as in text. 



direction. Interference between orbital and Zeeman coul- 
ing thus leads to a current in the negative y-dircction. 



.•(0) 



Quantitatively, we estimate from Eq. ( 10 ) jy, 

10 /iA/m and j^?'' ~ 1 nA/m for the current densities 
parallel and perpendicular to the plane of incidence, re- 
spectively, using a laser power of 1 W/mm^ as well as 
the parameters g = gs = 1 and vz = %^ ~ 29 m/s. For 
normal incidence the current vanishes, and for oblique 
incidence the dominant response is in the negative y'- 
direction since vz ~ 10-5 (cf. Fig.[l|. 

Results for more realistic dispersion. — When we in- 
clude deviations from the perfect Dirac cone in Eq. ^, 
there are additional hclicity-dependent and independent 
contributions to the photocurrent, as listed in Tabs. [l] 
and HD Without an external magnetic field, the leading 
correction to the helicity-dependent current is also in the 
x'-direction perpendicular to the plane of incidence. It is 
given by jII^I, (see Tab. |l| with 10"' 
and results from pure Zeeman coupling. There are no 
helicity-independent corrections to the photocurrent j^?' , 
so while the leading response in the a;'-(2/'-)direction is 
helicity dependent (independent), the overall leading re- 
sponse is helicity independent and parallel to the plane 
of incidence. 

Including an external magnetic field of strength 
B = 1 T, the leading helicity-dependent correction is 
•?hd with components parallel and perpendicular to 
the plane of incidence. The relative magnitude is 

ihd^Vii^'' ~ ^B/vz ~ 10"^ and results from interfer- 
ence between orbital and Zeeman coupling. The leading 
helicity-independent correction is j^^^^^ also with com- 
ponents parallel and perpendicular to the plane of inci- 
dence. This results solely from orbital coupling and has 

a relative magnitude of j^f^^Viy'*'' BA/vz ~ lO"'^. 
If the plane of incidence does not coincide with the yz- 
plane (cf. Fig. [T]), the dominant response both parallel 
and perpendicular to the plane of incidence is helicity 



independent with the relative magnitude of the perpen- 
dicular component x'/jhff'^^ ^z/BA ~ 10"'^. 

For normal incidence, i.e., 9 = if — 0, j^^'^ vanishes 
and there is no photocurrent, helicity-dependent or in- 
dependent, in the absence of an external magnetic field. 
Since j^?'' even vanishes when including B, the helicity- 
dependent photocurrent for normal incidence is always 
significantly smaller than for oblique incidence. The 
leading contributions are j^l^^J^ and jl^^J^ parallel to 



the magnetic field with j'hd k^VjI"^ 



10- 



10"'- Jm"^'^ is the only 



BA/vz 

and jH^'/j^ - XB^A/vl^ 

helicity-dependent contribution induced by pure orbital 
coupling. The leading helicity-independent contribution 
again is j'^f^^^ perpendicular to the magnetic field. For 
normal incidence including a magnetic field, the leading 
response in the si-direction is helicity dependent, while 
the y-direction response is helicity independent. 

The different contributions to the current as listed in 
Tabs. [I| and [II] depend strongly on the angle of incidence 
of the laser and thus the relative magnitude can change 
significantly with azimuthal angle tp of the vector poten- 
tial. The dependence of the current on ip tor 9 = 0.98 is 
plotted in Fig. [2] The response j^*^^ is independent of the 
polar angle but the corrections from B and A show strong 
angular dependence. The dominant current in the y'- 
direction is not affected by changes in the azimuthal an- 
gle but the dominant current in the x'-direction changes 
significantly as mentioned above. The large helicity- 
independent current jl^^i^^ vanishes for light incident in 
the yz-plane. 

Conclusion. — Motivated by recent experiments, we 
studied photocurrents in topological insulators. Focusing 
on the photocurrent response of the topological surface 
states, we find that the dominant photogalvanic current 
induced by obliquely incident circularly polarized light is 
helicity independent and in the plane of incidence of the 
light. This contribution is the result of an interference 
effect between the orbital and the Zeeman coupling of 
the light to the surface electrons. The helicity-dependent 
photocurrent is found to be very small. Although pure 
orbital coupling is the dominant excitation process, it 
does not induce a net photogalvanic charge current origi- 
nating from the surface states unless when including both 
band curvature and an in-plane magnetic field. Our re- 
sults suggest that an understanding of the experiments 
may require one to extend the theory to include the bulk 
states, the photon drag effect, or thermoelectric effects 
originating from inhomogeneous laser excitation. 
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Appendix A: Calculational Details 



When calculating the photocurrent density induced on the surface of a topological insulator given by (Eq. ^ in 
the main text), 



3 = 



^^^p,+|(p,+|iii|p,-)p5(2£;-M, 



(Al) 



one needs to be careful to include the perturbations on the Dirac spectrum in the interaction matrix element, as well 
as in the velocity and the ^-function. 

This section is organized as follows. First we will show how to calculate the integrand Vp_+\{p^+\H'_\p,—)\'^. 
Second, we will address how to treat the Delta function before calculating the entire integral in Eq. ( Al ). 



The calculation in detail of the integrand in Eq. ( Al ) proceeds as follows. Using projection operators onto the two 
bands defined as 



P.:Hp,±)(p,±|4±^, 
where is the full energy with all corrections to the perfect Dirac cone, the integrand can be expressed as 



(A2) 



rif-f 

| (p, +\H'_ |p, -) |2 = (p, Ip, +) (p, +\H'_ |p, -) (P, -I ^ |P, -) = tr 



P^H'P+H'_P- 



dH 
dp 



(A3) 



In the first step we used that the transition |p, +) — > |p, — ) {\p,—) \Pi+)) happens via emission (absorption) 
of a photon and is therefore mediated by the coupling involving (A_). The minus sign is included by taking 
the expectation value of the velocity operator with respect to the lower band. In the second step we used that the 
expectation value of an operator can be written as 



(p,±|0|p,±)=trIP±OI = tr 



i±^)o 

2 2E J 



(A4) 



In order to calculate Eq. ( A3 ) explicitly, we need to compute the interaction Hamiltonian given by (Eq. ([s]) in the 
main text) 



H' = —eVp^+ ■ A — gs^LBiy 'k A) ■ a. 



(A5) 



Since we would like to separate terms which create/annihilate a photon, i.e., write the interaction Hamiltonian in 
terms of A± , we write the velocities as 



dH dH dp+ dH dp^ 



w+ + V- 



and 



dH dH dp+ dH dp- 



dpx dp+ dpx (9p_ dpj: ^ dpy dp+ dpy dp^ dpy 

with p± ^ px i: ipy and dH/dp± = v± . The velocity operator then becomes 

dH dH ^ dH ^ , ,^ 

— = ^—x + ^—y^ [v+[x + ly) + v^{x - iy)\ . 
dp dpx dpy 



= «(w+ — W-). 



(A6) 



(A7) 



Here v± contains all corrections to the perfect Dirac cone. The vector potential for left circular polarization, k 
sin 9y — cos 6 z with e [0, 7r/2], and 93 = is given by 



A{t) — Ao[cos(fc • r — ijjt)x — sin(fe • r — Lot)(x x fc)] 

= -[A+e"''°'''(& - icos{9)y ~ism{0)z) + A^e'^''"^{x + icos{e)y + isin(6l)i) 



(A8) 



with A± — j4oe^*"*. From this we can read that 

1 



A, 



^{A++A_) and = — cos6'(A+ - A_). 



(A9) 
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With these expressions the orbital part of the interaction Hamihonian becomes 



dH OH 1 i 

-A,{t) + j—Ay{t) = v-)iA+ + A.) + -iv+~ v,)i~i cos 9){A+ - A.) 



dpx 



dpy 



-{A+ + C0S6') - cose*)] + A_ [u+(l - cos 6*) + w_ (1 + cos 6*)] }, (AlO) 



with H± being the part proportional to A± . In order to calculate the Zeeman term in the interaction using Eq. ( A8 ) 
we can calculate 

V X A = ^{^+[ - ikx {x~ icos{9)y - isin{9)z)]e^'''"^ + A^ [ik x {x + icos{9)y + isin{9)z)]e''^"^^ 
= ^ |^+ [ - x + i cos{9)y + i sin(6l)i] e"''"''" + A^[ - x - i coa{9)y - i sin(6')£] e*'"'''| 



-kA. 



(All) 



Note that for right circularly polarized light we would get the opposite sign in the last line of Eq. ( All[ ),i.e. , 
V X Arcp = ^^RCP- With Eq. (lAllI) the Zeeman coupling can be written as 



(V X A) • o- = [^^(crj, — i cos 9ay — i sin 6*0-2) + A^^a^ + i cos 9ay + i sin 0(7 ^ 



(A12) 



Inserting Eqs. (AlO) and (A12| into the interaction Hamiltonian given by Eq. ([3| in the main text, we get 
H' = - [u+(1 + cos6i) + w_(1 - cos 6*)] + A_ [w+(l - cos6l) + w_(l + cos6i)] 

— vz [A^ {ax — i cos 9ay — i sin 9a z ) + A^ (a^ + i cos Oay + i sin 9az)\ | , 



(A13) 



where vz 



is helicity dependent. The integrand can now be explicitly written as 



— t>p_+|(p, +\H'\p, — )p =— tr|P_A+ [w+(l + cos9) +v_{l — cos 9) — vzic^x — i cos 9ay — i sin 6'a-2)] 

P+A^ [v+{l — cos9) + v^{l + cos9) — vz{(Jx + icos9ay + i sin 6*0-2)] P_ 
[v+{x + iy) + v^{x- iy)] |, 



(AM) 



where the trace is understood to be over matrices only, not including the unit vectors. Eq. ( A14| is exact and contains 



all contributions from perturbations on the perfect Dirac spectrum arising from the interaction matrix element and 
the velocity. 



In addition to being careful with the definition of the velocity, we must also make sure the delta function expresses 
the perturbation on the Dirac spectrum. In order to Taylor-expand the argument of the (5-function, we use 



S{2E - huj) = AESiAE"^ - h^uj'^) = AE 



27r 



(A15) 



After expanding the entire integrand of the momentum integral in Eq. (All, i.e., Eq. (A14) and ( A15[ ), in the 
parameters A, A, and B to first order in A and A and second order in B, and performing the angular integral, we find 



j^'-^ Ida 



dp p 



S(0)(p)+aS(i)(p) + a2s(2)(p) 



(A16) 



where ^''■'(p) are functions of momentum and contain the parameters A, A, and B to the desired order. The exlicit 
expressions for the S(')(p) are very long without giving any insight and will not be presented here. The integral over 



a is simplified by writing a 



for the factors of a in the brackets. Integrating each summand separately. 



the derivatives can be pulled in front of the integrals. The integration over a in combination with the exponential 
function can now be resubstituted by a 5-function. Using S {4{pv p)"^ — — g^^'^^^p— 2t^^ ^^^^ ^ ~ (^)^j Eq. (A16) 
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Jp 


pre fact or 


X 


2/ 











A 




— sin(3(/)) 


cos(3(p) 


A 


^0 _ T T /I 

^vzXAcos6 


sin(3i/3) 


— cos(3(,3) 


B 


^vzXB^ COS0 


sin (/3 




ABi 


gBA 


— cos^ f> sin If 


3 cos^ 8 cos 


AB2 


fl4«A 


sin 


cos 1/3 


AB3 


ffilzAB^AcosS 


—9 sin ip + 5 sin(3(/3) 


cos tp[7 — 10 cos(2<(3)] 



TABLE III. P-polarization. vz = «z/uf ~ 10"^, A = \(fiujf/{v%) ~ 10"^, A = A{hu)f/(v'%) ~ 10"^, and B = gfisB/ihuj) ~ 
lO^^-B/T are dimensionless parameters and C is given in the main text. 





prefactor 


x' 


y' 





— ^vz s'mO 




1 


A 


^vzXcosB 


sin(3!^) 


— COs(3ly3) 


A 


^vzXAcose 


— sin(3(/3) 


COs(3(,5) 


Bi 


^vz\B^ cosO 






— COS ip 


B2 


§vlXBsm{20) 




— sm{2ip) 


cos(2(^) 


ABi 






— 3 sin ip 


cos 


AB2 


-^v^BAcos'^e 




sin tp 


cos ip 


AB3 


^vzXB'^AcosO 




[2 sin ip + 5 sin(3(/3)] 


9 cos 99 + 5 cos(3(/3) 


AS4 


i^i;|ABAsin(26») 




sm{2tp) 


— COs(2ly9) 



TABLE IV. S-polarized. Parameters as in Tab. [ml 



becomes 

= i!^^^ L/.(o) (m ^ ^ A.(i) (^yi^ _ ^.(2) ('^^ \ .A17) 

This can be easily evaluated leading to the results given in the main text. 



Appendix B: Linear Polarization 

For completeness we also calculated the response for P- and S-linearly polarized light. For k = sin By — cos 9 z and 
(j) = the vector potentials for P- and S-polarized light are given by 

Ap = Aq cos{k ■ r — aji)[cos Oy + sin6' (Bl) 

and 

As ~ Aq cos(fc • r — ujt) x. (B2) 

The calculation of the photocurrcnt proceeds as for circularly polarized light and the results are listed in Tabs. |III| 
and |IV| The resulting currents are, of course, helicity independent and the sum of the contributions from S- and P- 
polarized light add up to the helicity independent photocurrent induced by circular polarized light. As mentioned in 
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the main text, the overall largest contribution to the photocurrent, which is helicity independent and in the direction 
opposite to the direction of propagation of the incident light, can by induced by S-polarized light. 
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